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1.  Introduction 

The  airborne  Gravity  Gradiometer  Survey  System  (GGSS)  took  the  first  test  flight  in  the 
Texas/Oklahoma  area,  and  the  first  results  were  published  by  (Brzezowski  et  al.,  1988).  The  data 
released  by  the  Air  Force  Geophysics  Laboratory  (Jekeli,  circular  letter  5,  1989)  includes  the 
gradients  of  the  gravity  of  the  earth  measured  along  the  east-west  or  south-north  direction  in  19 
tracks.  The  GGSS  gradient  data  are  the  reduced  gravity  gradients.  The  gravity  gradients  were 
measured  by  the  GGSS  and  then  the  mean  values  of  the  measurements  of  each  track  were 
subtracted.  It  will  be  shown  later  that  these  reduced  gravity  gradients  are  the  second-order 
derivatives  of  the  disturbing  potential  plus  biases.  Therefore  it  is  important  to  remove  the  biases  in 
the  GGSS  gradient  data,  before  they  are  used  for  recovering  the  gravity  disturbance. 

The  GGSS  flew  along  the  east-west  (west-east)  or  south-north  (north-south)  direction  in  the 
test  area.  This  mission  supplies  the  crossover  discrepancies  of  the  gravity  gradients.  These  data 
are  useful  for  the  determination  of  the  biases  in  the  GGSS  gradient  data.  Like  the  determination  of 
the  bias  in  the  satellite  altimeter  data  for  each  arc,  the  relative  bias  in  the  GGSS  gradient  data  can  be 
determined  by  adding  some  constraints  to  the  crossover  discrepancy.  After  the  adjustment,  the 
biases  between  the  tracks  are  removed. 

The  recovery  of  the  first-order  derivatives  of  the  disturbing  potential  from  the  second-order 
derivatives  is  in  principle  an  integration  procedure.  In  this  procedure  the  absolute  bias  in  the 
GGSS  gradient  data  need  to  be  removed,  but  not  only  the  relative  bias.  The  absolute  bias  in  the 
GGSS  gradient  data  can  be  determined,  provided  tie  point  values  of  the  gravity  gradient  are  given 
at  the  flight  altitude.  Actually  no  such  tie  point  values  are  available,  so  the  absolute  bias  cannot  be 
determined  and  removed  from  the  GGSS  gradient  data.  Furthermore,  the  GGSS  data  varies  quite 
fast  in  the  area,  therefore  the  crossover  discrepancy  are  sensitive  to  the  positioning  error  of  the 
GGSS.  Based  on  the  above  consideration  we  would  say  that  the  crossover  discrepancy  adjustment 
is  not  suitable  for  the  removing  of  the  biases  in  the  GGSS  gradient  data.  But  this  adjustment  has 
another  application:  it  indicates  the  accuracy  of  the  measurement  of  the  GGSS  gradient  data.  This 
will  be  shown  later  in  the  paper. 

In  this  paper  we  use  linear  integration  to  transform  the  components  of  gravity  gradient  TXy, 
Tyz  into  the  vertical  derivative  of  the  disturbing  potential.  Then  the  tie  point  values  of  the  free-air 
anomaly  on  the  ground  are  combined  in  a  crossover  discrepancy  adjustment  for  the  determination 
of  the  bias  and  tilt  of  each  track.  Finally,  the  vertical  derivative  of  the  disturbing  potential  is 
determined  on  the  ground  by  using  the  fast  Fourier  transformation  (FFT)  and  least  squares 
collocation  (LSC).  The  results  are  compared  with  the  4  km  x  4  km  point  free-air  gravity  anomaly 
in  the  test  area  (Rapp  et  al.,  1988). 


2.  Variation  of  the  gradient  of  the  normal  gravity  in  a  local  area 

The  gravitational  potential  of  the  earth  W  is  usually  split  into  the  normal  and  disturbing 
potential  (Heiskanen  and  Moritz,  1967,  p.  82): 

W  =  U  +  T  (1) 

where  U  is  the  normal  potential  and  T  is  the  disturbing  potential. 

The  local  coordinate  system  xyz  (North-East-Down)  is  often  used  for  the  GGSS.  The  gravity 
gradient  Wjj  is  measured  in  the  xyz  coordinate  system  and  has  the  form: 

Wjj  =  U,j  +  T.j  i,  j  =  1, 2. 3  (2) 


1 


where  the  subscripts  i,  j  denote  the  derivatives  of  the  potential  and  the  numbers  1,  2  and  3 
corresponding  to  the  coordinates  x,  y,  z  respectively,  e.g.. 


Wj3  = 


2 

d  w 

dxdz 


=  WT 


(3) 


and  so  on. 


Because  the  GGSS  measures  only  the  gravity  gradient  Wy,  we  need  to  know  the  normal 
gravity  gradient  Uy  and  this  can  always  be  done  by  taking  account  of  the  reference  gravity  field. 

The  reference  gravity  field  is  chosen  as  a  normal  gravity  field  (ibid,  section  2-7): 

1  2  2  2 

U  =  V  +  i©  r  cos  w 

2  (4) 


where  V  is  the  normal  gravitational  potential;  to  is  the  angular  velocity  of  the  earth's  rotation;  r  is 
the  distance  between  the  computed  point  and  the  geocenter  of  the  earth  and  0  is  the  polar  distance 
of  the  computed  point  from  the  north  pole  (cf.  ibid.,  p.  18,  Figure  1-7). 

As  an  approximation  we  can  take  the  normal  potential  as  follows  (cf.  ibid,  eq.  (2-92)  or  (2- 
134a)): 


U  =  ^  -  J2  kM  cos  20  +  IJ  + 1  to2  r2  cos2  0 


where  k  is  the  gravitational  constant  and  M  the  mass  of  the  earth. 

In  the  spherical  coordinate  system  we  have  the  derivatives  of  the  normal 
order: 


(5) 

potential  in  the  first 


kM 

Ur  =  -  -~  +  3  J2kM  ^-l^-cos  20  +  ^-1+  co  rcos‘0 


a2 13 


Uo  = 


3  .  ..  a  1  2  : 
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sin  20 
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Ux  =  0,  (6) 

and  the  derivatives  of  the  normal  potential  in  the  second  order  are  given  by: 

U,r=2^-  12J2kM^cos20  +  i)  +  <o2cos20 
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V 


U*=- 


1  2  2 

2  J 2  kM  ~  +  co  rjsin20 


Uee- 


t  2  2 

3J2kM^--<o  r2 
r  / 


cos  20 


UrX  =  0 


Uex  =  0 

Uxx  =  o  (7) 

Usually,  the  airborne  flight  is  at  an  almost  constant  altitude  above  the  sea  level,  and  the  flight 
area  is  not  exceeded  by  a  few  degrees  of  latitude  and  longitude.  Thus,  the  second-order  derivatives 
of  the  normal  potential  can  be  written  in  a  linear  form: 

U„=uJ+<Ar  +  U^Ae 


Urt=U^+UrtAr  +  U°eA0 

Uee=  U  ee  +  U  Ar  +  UqqqAO 
UrX  =  0 


Uex  =  o 

Uxx  =  0  (8) 


where  Ar  =  r  -  ro,  A0  =  0  -  0o,  ro  and  0o  are  the  coordinates  of  the  inertial  point  which  can  be 
defined  as  any  point  at  the  beginning  of  the  GGSS  mission.  The  superscript  "0"  denotes  the 
values  of  the  derivatives  of  the  normal  potential  at  the  inertial  point. 

The  non-zero  third-order  derivatives  of  the  normal  potential  are  given  by: 


6kM 
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Um=-^  +  6OJ2kM^cos20+^ 
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r 

2  2 
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cos  20 
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a2  2  2 

Ueee=-  6J2kM  — -2co  r  sin  20 
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The  relationship  between  the  local  coordinate  system  xyz  (which  has  the  inertial  point  as  the 
origin)  and  the  spherical  coordinate  system  i0A  is  given  by 

I  x  =  r0(e  -  e0) 

<  y  =  r0cos  0  (A  -  A0) 

I  z  =  r0  -  r 

\  v  / 1  r\\ 


where  (ro,  0o,  Ao)  is  the  coordinate  of  the  inertial  point. 

Because  the  differences  AX  =  A  -  Ao  and  A0  =  0  -  0o  are  small  for  a  local  area,  we  have 
approximately 


cos  0  =  cos  (0O  +  A0)=»  cos  0O  +  A0  sin  0O 


If  we  choose  A0  <  3°  =  0.05  rad,  then  the  last  term  in  (11)  can  also  be  omitted,  and  eq.  (10) 
becomes: 


x  =  r0(o  -  0O) 
y  -  r0cos  0O(X  -  A0) 


z  =  r0  -  r 


From  (12)  we  have  immediately 

2  2 

a  a  a  a 


3r  dz 


’  2  2  ' 
dr  dz 


2  2 

_a__j__a_  a  _  i  a 

3e  =  r»^  882  =  rSax2' 
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3X  r»cose»3y  ax'  (r0coSe0)2ay2 


By  using  eqs.  (8),  (9)  and  (13);  we  get  the  variation  of  the  gradient  of  the  normal  gravity  in  a 
local  area: 
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^Uyy“  Uyy-  Uyy=  0 


AUxy=Uxy-U^y  =  0 


AUyz=Uyz- 1)^  =  0 

AUxx=Uxx-U;x=l(uee-Uee) 
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IT0  Z  ..o  y 

“ '  uiee"^T+  ueee“ 
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AUIZ=  U^z*  Ujj 


-U„-U 


o 
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(14) 


where  we  have  used  z  =  ro  -  r  =  -  Ar. 

Eq.  (14)  shows  that  the  gradients  of  the  normal  gravity  do  not  change  with  longitude  (y).  But 
they  vary  with  the  latitude  (x)  and  the  altitude  (z). 

Now  we  give  a  concrete  example  to  demonstrate  the  variation  of  the  gradients  of  the  normal 
gravity  in  a  local  area.  Taking  the  values  of  the  parameters  (Moritz,  1988,  p.  355): 

J2=  108  263  x  10'8 
a  =  6378137  m 


kM  =  3986005  x  108  m3  S'2 


to  =  72921 15  x  10-11  rad  S1,  (15) 

and  assuming  that  ro  =  a  +  1000  m,  0o  =  30°,  Xo  =  0,  then  we  have 
U^.=  -  1.44  x  10  3E6tvos/m 
U^=  -  2.13  x  10 1  Eotvos 
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U^e  =  •  3.39  x  107  Eotvds  .  m 


U0e9  =  2.13  x  1013E5tvos  .  m2 


Putting  eq.  (16)  into  eq.  (14),  we  obtain 


AU  =0 
yy 


AU  =0 

xy 


AUy2=  0 


AU  xx=  8.34  x  10'7  z  +  8.9 1  x  10'8  x 


AUXZ=  3.34  x  10'6z+ 8.34  x  10'7x 


AUzz=  1.44  x  10  3z  -  3.34  x  10  6 x 


From  eq.  (17)  we  can  see  that  the  gradients  of  the  normal  gravity  change  a  very  little  in  a  local  area. 
For  example,  let  z  =  300  meters,  x  =  300  km,  we  then  have 


AUxx  =  2.67x  10  Eotvos 


AUXZ=  2.50x10  Eotvos 


AUn=  -  0.57  Eotvos 


From  eq.  (18)  we  can  see  that  only  AUzz  has  a  little  change  in  a  local  area,  and  other  gradients  of 
the  normal  gravity  can  be  considered  as  constant. 

Now  we  discuss  some  practical  aspects.  The  test  flight  of  a  GGSS  was  taken  in  a  local  area. 
If  the  gradient  data  are  collected  for  each  track  and  the  mean  value  is  subtracted  from  these  data, 
then  we  have  the  reduced  gradient  data  as: 

T?j=Wij.Wii  (19) 

where  the  overbar  denotes  themean  values  of  the  gradient  of  the  gravity  for  a  track,  Wjj  is  the 
observed  gravity  gradient  and  Tjj  is  the  reduced  gradient  of  the  gravity. 

Based  on  the  above  discussion  we  have  approximately 


Uij-Uj, 
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then  we  obtain 


Tij =  Wjj  -  (Ujj  +  TSJ 

•Wu-tr.j-T.j 

=  Tij "  Tij  (2i) 

Eq.  (2D  indicates  that  the  reduced  gradient  data  is  different  from  the  real  second-order 
derivaties  of  the  disturbing  potential  Tij  by  a  constant  Tjj,  Normally,  the  global  average  of  the 
disturbing  potential  and  its  derivatives  should  approach  zero,  but  the  mean  value  of  the  Tjj  of  each 
track  will  not  always  be  equal  to  zero.  This  is  one  of  the  causes  of  the  bias  in  the  GGSS  gradient 
data. 


3.  Relationship  Between  First-order  and  Second-order  Derivatives  of  the 
Disturbing  Potential 

The  relationship  between  the  first-order  and  second-order  derivatives  of  the  disturbing 
potential  is  given  by  (cf.  Moritz,  1971,  p.  13): 

Ti  =  T°  +  f  Pfrxidx  +  Tyidy  +  Tzidz)  i=  1,2,3 

p°  (22) 

Here  T^  is  the  initial  value  of  Ti  at  the  point  Po  and  Tj  refers  to  a  current  point  P.  Obviously,  the 
constant  bias  in  the  gravity  gradient  data  plays  an  important  role.  If  we  have  a  bias  c  in  the  data 
Txi,  then  from  eq.  (22)  we  have 


c . 


x 


(23) 


where  x  is  the  change  of  the  x-coordinate  of  the  curve  Po  P.  Eq.  (23)  means  that  a  bias  in  the 
gravity  gradient  data  becomes  a  tilt  in  the  gravity  disturbance.  Because  the  bias  is  different  for  each 
component  of  the  gravity  gradient  data  measured  on  each  track,  the  gravity  disturbances  obtained 
by  using  the  linear  integration  have  different  tilt.  Before  we  put  the  gravity  disturbances  on  each 
track  together  to  get  the  gravity  disturbances  on  the  ground,  the  tilt  of  each  track  has  to  be 
removed. 


4.  Crossover  Discrepancy  Adjustment  for  Processing  Aerial  Gravity  Gradient 
Data 


A  crossover  discrepancy  adjustment  has  been  used  for  processing  satellite  altimeter  data.  At 
first  glance  it  seems  that  the  crossover  discrepancy  adjustment  can  be  directly  used  to  the  gravity 
gradient  data  for  removing  the  bias  in  the  data.  As  in  processing  the  satellite  altimeter  data,  if  only 
the  crossover  discrepancy  and  some  constraints  are  used,  the  relative  bias  in  data  can  be 
determined.  The  relative  bias  determination  is  useless  for  gravity  gradient  data.  The  gravity 
gradient  has  six  components.  If  only  the  relative  bias  is  removed  from  each  component  of  the 
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gravity  gradient  for  each  track,  there  are  still  different  biases  in  components  of  the  gravity  gradient 
data.  The  reduced  data  cannot  be  combined  together  to  recover  the  gravity  disturbance  on  the 
ground.  Of  course  the  absolute  bias  in  the  gradient  data  can  be  determined,  provided  other  data  can 
be  used,  e.g.,  the  tie  point  values  of  the  gravity  gradient  at  the  flight  aldtutde,  but  this  kind  of  data 
does  not  exist  at  present.  Therefore  the  crossover  discrepancy  adjustment  is  not  practical  for  the 
elimination  the  bias  in  the  gradient  data. 

The  bias  in  the  gravity  gradient  data  is  a  serious  problem  for  recovering  the  gravity 
disturbance.  In  order  to  solve  this  problem,  we  transform  the  second-order  derivatives  of  the 
disturbing  potential  into  the  first  order.  The  absolute  bias  and  tilt  in  the  first-order  derivatives  of 
the  disturbing  potential  can  be  removed  by  using  the  crossover  discrepancy  adjustment  and  the  tie 
point  values  (gravity  anomalies,  deflections  of  the  vertical,  etc.). 

In  the  following  we  consider  using  two  components  of  second-order  derivatives  of  the 
disturbing  potential,  Txz  and  Tyz,  to  recover  the  vertical  component  of  the  gravity  disturbance. 

The  test  flight  of  the  GGSS  were  taken  along  south-north  (x)  and  west-east  (y)  directions. 
By  using  eq.  (22)  we  have 


Tz(x-yJ  =  Tz(xo-  yj-*-  /  Txz(x'*yj)dx' 

*0 

Tz|xi-  y)  =  Tz(Xi,  yo )  +  /y  T  (x;,  y)dy' 

y0 


(24) 


(25) 


where  the  subscripts  i,  j  denote  the  tracks  along  west-east  or  south-nonh  direction,  respectively. 

Assume  that  there  are  random  errors  and  a  constant  bias  in  the  gravity  gradient  data.  Then 
eqs.  (24)  and  (25)  become: 


Tz(x,yJ  +  V(x,yJ  =  Tz(x0,  yj  +  f  Txzfx',  yjdx' +  b(yj(x  -  Xo) 

x°  (26) 

Tz(xi.  y}+  V'K  y)  =  Tz(xit  yo )+  f  Tyz(xi,  yjdy'  +  bjxjjy-yo) 

y°  (27) 

where  V  and  V'  are  the  errors  in  the  Tz  caused  by  random  errors  in  the  gravity  gradient  data,  b  (yj) 
and  b  (xj)  are  the  bias  in  the  gradient  data. 


Denote 


-  f  Txz(x.yj)dx, 
*/*0 

f  Tyz(xit  y)dy. 
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bi  =  b  (xj),  bj  =  b  (yj) 

ai  =  Tz  (xj,  y0),  aj  =  T2  (xo,  yj), 

then  eqs.  (26)  and  (27)  have  the  values  at  the  point  (xj,  yj): 

ifj  +  V'fxj,  yj  =  ar  TzjXj,  yj  +  b;  (yj  -  yj 
lij +  v (xi*  yj  =  aj  *  Tz(xi«  yj  +  bj(xi  -  xo) 


(28) 


(29) 

(30) 


Eqs.  (29)  and  (30)  can  be  considered  as  observation  equations.  The  "observed"  quantities  l*j 
and  lfj,  can  be  obtained  directly  by  using  the  numerical  integration  of  the  components  of  the  GGSS 
gradient  data  and  Tyz  Subtract  (30)  from  (29),  to  obtain  the  crossover  discrepancy: 

dij  +  V,j  =  lfj -  lfj  =  -  aj  +  bijyj  -  yo)  -  bj (xs -  xq)  (3 1} 


with  Vy  =  V’  (xj,  yj)  -  V  (xj,  yj),  dy  =  lfj  -  lfj. 

Eq.  (31)  is  the  basic  equation  for  the  crossover  discrepancy  adjustment.  The  crossover 
discrepancy  dy  is  given,  the  parameters  aj,  aj,  b;,  and  bj  will  be  determined.  Like  the  crossover 
discrepancy  adjustment  for  altimeter  data,  eq.  (31)  describes  an  indeterminate  problem,  because 
there  is  no  track  fixed.  In  order  to  remove  the  indeterminancy  the  tie  point  values  are  introduced. 

Now  we  consider  how  many  tie  point  values  are  needed  for  the  crossover  discrepancy 
adjustment.  If  there  are  only  biases  in  eqs.  (29)  and  (30),  i.e.,  bj  and  bj  are  equal  to  zero,  in 
principle  one  tie  point  value  of  the  Tz  is  needed  for  the  determination  of  absolute  bias  of  each  track. 
If  the  bias  and  tilt  are  determined  for  each  track,  three  tie  point  values  whose  locations  are  not  on 
the  same  line  are  needed  at  least.  This  is  the  same  as  the  theorem  of  the  determination  of  a  plane: 
Three  points  not  on  the  same  line  determine  a  plane. 

We  have  assumed  that  the  tie  point  values  are  errorless.  In  practice  it  is  not  the  case. 
Therefore  using  more  tie  point  values  is  more  reliable  than  using  merely  necessary  tie  point  values. 

If  we  put  the  tie  point  values  at  the  beginning  and  end  of  the  tracks  along  south-north 
direction,  by  using  eq.  (29)  we  get: 

di  +  b|(yk-yo) 

i  =  1,  2,  3...;  k  =  1,  2  (32) 


with 


+  Ifk •  V?  =  v'(*i.  yj  •  Vj(Xj,  yj  (33) 

where  (Tz)j  is  the  tie  point  value,  Vj  is  the  measurement  error,  yk  is  the  y-coordinate  of  the  tie 
point. 
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The  parameters  a;,  aj,  bj,  bj  can  be  determined  by  using  least  squares  adjustment  by  using  eq. 
(31)  and  (32)  together. 

Eqs.  (31)  and  (32)  can  be  written  in  matrix  form: 

D  +  V  =  AX  (34) 

where 


A  least  squares  solution  is  given  by 


T 

A  PD 


(35) 


where  At  is  the  transpose  of  the  matrix  A  and  P  is  the  weight  matrix.  P  is  chosen  as  a  diagonal 
matrix  which  has  the  elements  of  the  reciprocal  of  the  standard  error  of  the  discrepancy  djj  and  d* . 


5.  Downward  Continuation  of  the  Gravity  Disturbance  from  the  Flight  Altitude  to 
the  Ground  by  Using  FFT  and  Least  Squares  Collocation 

Suppose  that  the  bias  and  tilt  of  die  vertical  derivatives  of  the  disturbing  potential  are  removed 
by  using  the  procedure  from  section  4,  then  the  reduced  vertical  component  of  the  gravity 
disturbance  has  to  be  downward  continued  to  the  ground.  For  the  downward  continuation  we  can 
use  the  fast  Fourier  transformation  (FFT)  or  least  squares  collocation  (LSC). 

The  basic  formula  for  the  upward  continuation  problem  is  the  Poisson's  integral.  In  a  planar 
approximation  it  is  (cf.  Heiskanen  and  Moritz,  1967,  p.  240): 
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where  H  is  the  flight  altitude  of  the  GGSS,  d  is  the  distance  between  the  current  point  on  the 
ground  and  the  computed  point  at  the  flight  altitude.  Tz  and  Yz  are  the  vertical  components  of  the 
gravity  disturbance  at  the  flight  altitude  and  on  the  ground,  respectively. 


Applying  the  Fourier  transformation  to  eq.  (36),  we  get 
F  {Tz)  =  e’2n<oH  f{T;} 

where 

2  T 

U  +  V  ,  U,  V 

are  the  frequency  variables;  F  {  }  denotes  the  Fourier  transformation. 
From  eq.  (37)  we  have 

F{T;)  =  e2"“HF|TI| 


(37) 


(37*) 


Eq.  (37*)  describes  the  downward  continuation  of  the  Tz  to  the  ground,  and  it  is  used  often  for 
determining  the  vertical  component  of  the  gravity  disturbance  on  the  ground.  In  order  to  reduce  the 
effect  of  the  high  frequencies  of  the  error  in  downward  continuation,  some  kinds  of  smoothing  or 
filtering  should  be  used. 

Least  squares  collocation  is  used  in  two  aspects.  First,  there  are  only  a  few  tie  point  values 
given  on  the  ground.  They  are  often  too  sparse  to  use  eq.  (36)  to  be  upward  continued  to  the  flight 
altitude.  This  can  be  done  only  by  using  least  squares  collocation. 

Second,  least  squares  collocation  can  also  be  used  for  the  downward  continuation  of  the 
derivatives  of  the  disturbing  potential  from  flight  altitude  to  the  ground.  The  basic  formula  for  least 
squares  collocation  is  given  by  (Moritz,  1980,  p.  102): 


-l 


s  =  Csl(Ctl  +  Cnn)  t 


(38) 


where  t  is  the  observation,  Cst  is  the  covariance  between  the  computed  quantity  s  and  the 
observation  t,  Cu  is  the  covariance  between  the  observations,  Cnn  is  the  covariance  of  the 
measurement  error. 

It  is  reasonable  to  choose  the  covariance  function  for  the  vertical  component  of  the  gravity 
disturbance  as  (cf.  Moritz,  1971,  p.  36;  and  1980,  p.  179): 


C(P,Q)  = 


C0  b2(z  +  z'+  b) 

n3/2 

2  2 
p  +  (z  +  z'+b) 


(39) 


where  b  is  a  constant,  P,  Q  have  the  coordinates  (x,  y,  z)  and  (x',  y',  z'),  p  is  the  distance  between 
the  points  (x,  y,  0)  and  (x',  y\  0),  the  variance  Co  is  the  value  of  the  covariance  function  C  (p,  Q) 
for  P  =  Q. 
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The  airborne  flights  are  at  a  constant  altitude,  so  we  have  z  =  0,  z'  =  constant.  Eq.  (39) 
becomes: 


C(P,Q)  = 


C0b2(z'  +  b) 

,2  +  (z'  +  b)2f 


The  covariance  function  will  be  used  for  the  upward  and  downward  continuation  of  a  harmonic 
function  above  the  sea  level. 

6.  Practical  Consideration  and  Computation  Procedure  Description 

The  GGSS  data  were  measured  along  the  tracks  oriented  west-east  and  south-north.  The 
mathematical  study  of  this  procedure  is  given  in  section  2  and  the  conclusion  is  that  the  data  can  be 
considered  as  the  second-order  derivatives  of  the  disturbing  potential  plus  a  bias.  Therefore,  the 
most  important  problem  in  processing  the  GGSS  gradient  data  is  the  bias  problem. 

By  using  the  crossover  discrepancy  adjustment  the  relative  bias  in  the  components  of  the 
gravity  gradient  can  be  determined.  If  the  absolute  bias  is  to  be  determined,  at  least  one  tie  point  of 
the  gravity  gradients  at  the  flight  altitude  is  needed. 

If  there  are  different  biases  in  the  components  of  the  second-order  derivatives  of  the  disturbing 
potential,  we  cannot  combine  all  six  components  together  to  recover  the  gravity  disturbance 
because  the  biases  make  the  gradient  data  inconsistent  for  this  recovery.  On  the  other  hand,  the 
crossover  discrepancy  indicates  the  accuracy  of  the  data.  Usually  the  measurement  error  (e.g., 
positioning  error,  instrument  error)  can  be  split  into  systematic  and  random  parts.  In  the  same 
manner  the  errors  in  the  gradient  data  can  be  split  into  systematic  error  which  consists  of  the 
systematic  measurement  error  and  the  bias  for  each  track  (cf.  section  2),  and  random  error.  After 
the  crossover  discrepancy  adjustment,  the  bias  and  the  systematic  measurement  error  are  reduced, 
then  the  new  crossover  discrepancy  dij  is  equal  to  the  difference  of  the  random  pan  of  the 
measurement  error  in  the  gradient  data: 

NS  EW 

d-  •  =  p. .  —  f.  . 

»j  fc'j  (4n 


where  6jj ,  are  the  random  part  of  the  measurement  errors  of  the  gravity  gradients  along  tracks 
oriented  to  north-south  and  east-west  directions,  the  subscripts  i,  j  denote  the  values  and  are 
located  at  the  crossover  of  the  tracks  i  and  j. 

Assume  that  the  random  measurement  errors  e?jS  and  are  unconelated,  then  the  RMS 
value  of  the  discrepancy  djj  is  equal  to 


where  a?jS  and  are  the  standard  random  measurement  errors  of  the  tracks  along  north-south 
and  east- west  directions;  C?j  is  the  RMS  value  of  the  crossover  discrepancy  between  tracks  i  and  j. 


Therefore,  the  RMS  value  of  the  crossover  discrepancy  indicates  the  measurement  accuracy  of  the 
gradient  data. 

In  Table  1  we  give  the  statistics  of  the  crossover  discrepancy  for  the  GGSS  gradient  data 
before  and  after  the  crossover  discrepancy  adjustment.  This  table  shows  that  the  accuracy  of  the 
GGSS  gradient  data  is  around  10  Eotvos. 


Table  1 

Statistics  of  the  Crossover  Discrepancy  of  the  GGSS  Gradient  Data 

Unit  E6tv5s 


Txx 

~'T— 

'  T~ 

. ~Tvz 

— T^T 

mean 

-3.8 

4.4 

-2.4 

-3.6 

7.4 

before 

max 

37.5 

35.2 

28.3 

38.8 

38.0 

58.6 

adjust. 

min 

-36.9 

-32.1 

-36.1 

-48.3 

-55.6 

-40.5 

RMS 

14.4 

13.1 

12.1 

16.6 

17.9 

19.4 

STD 

13.9 

12.3 

11.9 

16.2 

17.9 

17.9 

mean 

-0.0 

-0.0 

-0.0 

-0.0 

0.0 

0.0 

after 

max 

34.6 

41.7 

27.9 

30.6 

54.9 

32.9 

adjust. 

min 

-28.6 

-29.5 

-26.1 

-32.0 

37.7 

-32.4 

RMS 

12.2 

11.2 

9.7 

12.6 

13.9 

14.2 

STD 

12.2 

11.2 

9.7 

12.6 

13.9 

14.2 

Usually  the  gravity  anomaly  rather  than  the  gravity  disturbance  is  given  in  an  area.  The 
relationship  between  the  free-air  anomaly  and  the  vertical  component  of  the  gravity  disturbance  is 
given  by  (Heiskanen  and  Mortiz,  1967,  p.  89): 


2T 

Ag  =  T*-TT 


(43) 


Notice  that  the  axis  z  is  in  the  opposite  direction  of  the  radius  vector  r.  As  first  approximation  we 
ignore  the  last  term  in  eq.  (43)  and  do  not  distinguish  the  free-air  anomaly  and  the  vertical 
component  of  the  gravity  disturbance. 


Tn  the  test  area  (Texas/Oklahoma)  4  km  x  4  km  point  free-air  anomaly  data  are  available. 
Therefore  we  transfonmed  the  components  of  the  gradient  data  Txz  and  Tyz  into  the  vertical 
component  of  the  gravity  disturbance  Tz  by  using  the  linear  integration.  Then  the  tie  point  values 
of  the  free-air  anomaly  were  used  to  determine  the  bias  and  tilt  in  the  Tz  for  every  track.  The 
computation  procedure  is  described  as  follows. 


a.  Taking  linear  integration  of  the  components  Tyz  and  Txz  along  south-north  (x)  and  west-east 
(y)  directions,  the  formulas  (24)  and  (25)  were  used; 


b.  The  4  km  x  4  km  point  free-air  anomaly  are  given  in  the  test  area.  Let  the  tie  points  have  the 
coordinates  (xj,  yj,  zo),  i  =  1,  2,...  20,  and  zq  =  1  km.  Here  xj,  yj  are  the  coordinates  of  the  tie 
points  located  at  the  beginning  and  end  of  the  tracks  along  the  south-north  direction  in  the  square 
area  (cf.  Figure  1).  The  4  km  x  4  km  point  free-air  gravity  anomaly  was  interpolated  at  the  points 
which  have  the  coordinates  (xj,  yj)  by  using  the  bicubic  spline  function,  then  was  upward 
continued  to  the  flight  altitude  by  using  least  squares  collocation.  The  covariance  function  used  is 
given  by  eq.  (39). 
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c.  Using  the  tie  point  values  of  the  free-air  gravity  anomaly  and  the  crossover  discrepancy  to 
determine  the  bias  and  tilt  for  each  track;  the  least  squares  adjustment  (eq.  (35))  was  used; 

d.  Using  the  bicubic  spline  function  to  interpolate  the  adjusted  Tz  at  the  flight  altitude  into  a  2  km 
x  2  km  grid  interval,  then  fast  Fourier  transformation  was  used.  As  a  comparison  least  squares 
collocation  was  used. 


7.  Results 

The  abridged  airborne  GGSS  tracks  are  drawn  in  Figure  1.  The  data  used  for  the  analysis  are 
constrained  in  the  square  area  circled  by  thick  lines.  There  are  19  tracks  in  which  9  tracks  are 
along  the  east-west  direction  and  10  tracks  are  along  the  south-north  direction.  After  the  linear 
integration  of  the  Txz  and  Tyz  along  the  south-north  (x)  and  west-east  (y)  directions,  1* j,  l?j  (cf. 
eqs.  (24)  -  (28))  are  given  along  each  track.  By  combining  the  tie  point  values  of  the  free-air 
gravity  anomaly  and  the  crossover  discrepancy  the  absolute  bias  and  tilt  were  determined  for  every 
track. 

In  the  least  squares  adjustment  solution,  eq.  (35)  the  weight  matrix  P  has  the  elements  (102 
mgal2)*1  for  the  crossover  discrepancies  and  (22  mgal2)*1  for  the  tie  point  values.  The  statistics  of 
the  crossover  discrepancy  are  given  in  Table  2. 


Table  2. 

Statistics  of  Crossover  Discrepancy  of  the  Tz 
Unit  mgal 


max 

min 

mean 

RMS 

STD 

Before  adj. 
After  adi. 

82.56 

16.62 

-403.60 

-17.55 

-109.13 

-0.00 

157.39 

7.75 

113.39 

7.75 

From  Table  2  we  can  see  that  the  crossover  discrepancy  is  reduced  enormously  after  the 
adjustment.  Based  on  eq.  (42)  one  can  assess  the  accuracy  of  the  Tz  as  about  5-8  mgal. 

After  the  crossover  discrepancy  adjustment,  Tz  is  given  at  the  flight  altitude.  It  can  be 
analytically  downward  continued  to  the  ground  by  using  the  fast  Fourier  transformation  or  least 
squares  collocation.  As  a  comparison  both  methods  were  used  and  the  results  are  given  in  Table  3 
and  the  following  figures. 

For  the  FFT  computation  a  2  km  grid  interval  was  used,  then  the  results  were  averaged  into  4 
km  x  4  km  mean  block  values.  For  the  least  squares  collocation  computations  the  4  km  grid 
interval  was  used,  and  in  the  LSC  solution  230  data  points  were  used  for  the  W-E  tracks  and  260 
data  points  for  the  S-N  tracks. 

Table  3  shows  the  statistics  of  the  differences  between  the  free-air  anomaly  Ag  and  the  vertical 
derivatives  of  the  disturbing  potential  Tz  which  is  recovered  by  processing  the  GGSS  gradient 
data.  This  comparison  may  be  somehow  questionable,  because  the  firee-air  anoamaly  and  the  Tz 
arc  not  the  same  quantities.  Even  so,  the  difference  between  them  (the  last  term  in  eq.  (43))  is 
small  in  comparison  with  the  free-air  anomaly  and  Table  3  can  still  give  us  a  general  idea  about  the 
quality  of  the  recovering  of  the  Tz  from  the  GGSS  data. 
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Tabic  3 


Statistics  of  the  Differences  Between  Ag  and  Tz  in  Test  Area 

Unit  mgal 


wnmmm 

max 

■32m 

mean 

RMS 

STD 

Ag  -  Tz  (W-E) 

FFT 

LSC 

-17.54 

-18.93 

2.47 

0.85 

8.60 

8.04 

8.26 

7.99 

Ag  -  Tz  (S-N) 

FFT 

LSC 

27.76 

23.71 

-20.19 

-19.65 

2.27 

2.29 

8.75 

7.96 

8.45 

7.62 

Ag  -  Tz  (all) 

FFT 

LSC 

24.14 

14.73 

-20.46 

-16.93 

-1.77 

1.49 

7.70 

5.94 

7.49 

5.75 

In  Table  3  the  Tz  (W-E)  and  Tz  (S-N)  are  the  Tz  obtained  by  using  the  data  on  tracks  along 
west  -east  and  south-north  directions,  respectively;  Tz  (all)  is  the  Tz  obtained  by  combining  all  Tz 
data  on  the  flight  altitude  together. 

Table  3  shows  that  both  methods  FFT  and  LSC  are  given  almost  the  same  results.  The  results 
are  better  if  least  squares  collocation  is  used.  The  reason  is  that  least  squares  collocation  gives  a 
smoother  result  than  the  fast  Fourier  transformation  does.  This  will  be  clearly  shown  in  the 
following  Figures. 

Figure  1  shows  the  abridged  GGSS  tracks.  The  data  used  for  the  analysis  are  confined  inside 
the  square  area. 

Figure  2  is  the  contour  map  of  the  free-air  anomaly  based  on  the  4  km  x  4  km  grid  in  the  test 
area.  The  free-air  anomaly  varies  significant  along  direction  south-west  to  north-east. 

Figure  3  and  4  shows  the  Tz  recovered  by  using  the  Tz  data  on  the  tracks  along  west-east 
directions.  Figure  3  is  the  FFT  solution  and  4  is  the  LSC  solution.  Obviously,  the  LSC  solution 
is  smoother  than  FFT  solution. 

Figure  5  and  6  are  the  contour  maps  of  the  Tz  obtained  by  using  the  Tz  data  on  the  tracks  along 
south-north  direction. 

Figure  7  and  8  show  the  Tz  recovered  by  using  the  Tz  data  on  the  tracks  along  west-east  and 
south-north  directions.  Obviously,  they  are  the  best  results. 
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Figure  2.  Contour  Map  of  the  Free-air  Anomaly  in  Test  Area 
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Figure  3.  Contour  Map  of  the  Tz  (W-E)  FFT  Solution 
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Figure  4.  Contour  Map  of  the  Tz  (W-E)  LSC  Solution 
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Figure  5.  Contour  Map  of  the  Tz  (S*N)  FFT  Solution 
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Figure  6.  Contour  Map  of  the  Tz  (S-N)  LSC  Solution 
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Figure  7.  Contour  Map  of  the  Tz  (all)  FFT  Solution 
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Figure  8.  Contour  Map  of  the  Tz  (all)  LSC  Solution 
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8.  Conclusions  and  Recommendations 


Because  there  are  biases  in  the  GGSS  gradient  data  for  different  components  of  the  gravity 
gradient  and  different  tracks,  this  data  cannot  be  directly  used  for  recovering  the  gravity 
disturbance  on  the  ground.  In  order  to  remove  the  bias  in  data,  the  crossover  discrepancy 
adjustment  and  constraint  equations  can  be  used.  By  using  the  crossover  discrepancy  adjustment 
and  constrain  conditions,  the  relative  bias  in  the  data  can  be  determined  and  removed.  Even  so  the 
gradient  data  are  difficult  to  be  combined  together  to  recover  the  gravity  disturbance,  because  the 
different  biases  exist  in  the  components  of  the  gravity  gradient  after  the  crossover  discrepancy 
adjustment  and  they  make  the  data  inconsistent  for  the  recovery  of  the  gravity  disturbance. 

In  this  study  we  used  the  linear  integration  to  transform  the  second-order  derivatives  of  the 
disturbing  potential  Txz  and  Tyz  into  the  vertical  derivative  of  the  disturbing  potential  Tz,  then  the 
crossover  discrepancy  and  the  tie  point  values  of  the  free-air  anomaly  were  used  to  remove  the  bias 
and  tilt  in  the  Tz.  The  best  result  agrees  with  the  free-air  anomaly  in  the  test  area  at  the  6  mgal  level 
(RMS  difference).  In  comparison  with  the  accuracy  of  the  GGSS  gradient  data  (the  crossover 
discrepancy  indicates  that  the  GGSS  gradient  data  could  have  the  accuracy  about  10  Eotvds),  the 
accuracy  of  the  Tz  is  reasonable. 

In  the  computations,  the  FFT  and  LSC  method  were  used;  both  gave  almost  the  same  results. 
The  LSC  solution  is  smoother  than  the  FFT  solution  and  agrees  a  little  better  with  the  free-air 
anomaly  in  the  test  area  (cf.  Table  3). 

The  weakness  of  the  computation  procedure  is  that  only  two  components  of  the  GGSS 
gradient  data  were  used.  In  the  crossover  discrepancy  adjustment  for  the  Tz  there  were  20  tie  point 
values  used. 

In  principle,  one  tie  point  (one  tie  point  has  six  components  of  the  gradient  data)  is  needed  to 
determine  the  absolute  bias  in  the  gradient  data.  If  such  tie  point  values  are  available,  the  absolute 
bias  in  the  gradient  data  can  be  determined  and  removed.  The  reduced  gradient  data  can  then  be 
combined  together  to  recover  the  three  components  of  the  gravity  disturbance.  Normally,  there  are 
measurement  errors  in  the  tie  point  values.  More  tie  point  values  can  also  improve  the  results  of 
the  crossover  discrepancy  adjustment  by  using  least  squares  adjustment.  In  the  future,  if  we  get 
the  tie  point  values  of  the  gravity  gradient,  then  the  gradient  data  can  be  used  together  to  get  better 
recovery  of  the  gravity  disturbance. 

If  the  second-order  derivatives  of  the  disturbing  potential  are  transformed  into  the  first  order, 
then  the  bias  and  tilt  have  to  be  determined  for  each  track.  In  principle,  three  tie  points  (every  tie 
point  has  the  Tx,  Ty  and  Tz  three  values)  are  needed  to  remove  the  bias  and  tilt  in  the  gravity 
disturbance  data.  If  the  accuracy  of  the  gradient  data  is  poor,  as  is  the  case  in  our  computations, 
more  tie  point  values  are  needed  to  improve  the  computation  accuracy. 

By  the  way,  the  GGSS  is  still  in  development.  The  accuracy  of  the  measurements  will  be 
improved.  One  can  hope  in  the  future  that  one  can  get  accurate  gravity  disturbance  on  the  ground 
by  processing  such  kinds  of  data.  If  the  GGSS  measurements  are  accurate  enough,  only  a  few  tie 
point  values  are  needed  and  this  makes  the  GGSS  more  efficient  and  useful. 


9.  References 

Brzezowski,  S.,  D.  Gleason,  J.  Goldstein,  W.  Heller,  C.  Jekeli,  and  J.  White,  Synopsis  of  Earth 
Field  Test  Results  from  the  Gravity  Gradiometer  Survey  System,  In:  Chapman  Conference 
on  Progress  in  the  Determination  of  the  Earth's  Gravity  Field,  Sept.  1988. 


24 


Heiskanen,  W.A.,  and  H.  Moritz,  Physical  Geodesy,  W.H.  Freeman,  San  Francisco,  1967. 

Jekeli,  C.,  Terrestrial  and  Airborne  Gravity  Gradiometry,  Circular  Letter  5,  1989,  I.A.G.  SSG 
3.111,  Circular  Letter  5,  1989. 

Moritz,  H.,  Kinematical  Geodesy  II,  Report  No.  165,  Department  of  Geodetic  Science  and 
Surveying,  Ohio  State  University,  Columbus,  1971. 

Moritz,  H.,  Advanced  Physical  Geodesy,  Herbert  Wichmann  Verlag,  Karlsruhe,  1980. 

Moritz,  H.,  Geodetic  Reference  System  1980,  Bulletin  Geodesique,  Vol.,  62,  No.  3,  1988. 

Rapp,  R.H.,  and  S.  Zhao,  The  4  km  x  4  km  Free-Air  Anomaly  File  for  the  Continuous  United 
States,  Internal  Report,  Dept,  of  Geodetic  Science  and  Surveying,  The  Ohio  State  University, 
1988. 


25 


